Decoherence in a Josephson junction qubit 
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The zero-voltage state of a Josephson junction biased with constant current consists of a 
set of metastable quantum energy levels. We probe the spacings of these levels by using microwave 
spectroscopy to enhance the escape rate to the voltage state. The widths of the resonances give a 
measurement of the coherence time of the two metastable states involved in the transitions. We 
observe a decoherence time shorter than that expected from dissipation alone in resonantly isolated 
20x5 (/xm) 2 area Al/AlOx/Al junctions at 60 mK. The data is well fit by a model that includes 
the dephasing effects of both low-frequency current noise and the escape rate to the voltage state. 
We discuss implications for quantum computation using current-biased Josephson junction qubits, 
including limits on the minimum number of levels needed in the well to obtain an acceptable error 
limit per gate. 

PACS numbers: 03.65.Yz,03.67.Lx,85.25.C,78.70.Gq 



Research in the 1980s definitively showed that 
the phase difference across a single current-biased 
Josephson junction can behave quantum-mechanically 
|l|,|2j. The recent proposal that an isolated current-biased 
Josephson junction could serve as a qubit Q in a quan- 
tum computer has preceeded a resurgence of interest in 
this simple system 0, 0,0,0 • 

Designing a quantum computer based on isolated 
Josephson junctions raises many issues. Isolation of the 
junction from its bias leads must be achieved by control- 
ling the high frequency electromagnetic environment that 
the junction couples to 2]. At the very least, this iso- 
lation must be effective around the resonant frequency 
of the junction. In addition, at lower frequencies, cur- 
rent noise will tend to cause decoherence in the junction 
state. Also, during typical gate operations the junction 
must operate in a strongly anharmonic regime that can 
be reached by applying a large bias current through the 
junction. In this high bias regime however, there is an 
increased escape rate from the upper qubit state. In this 
Letter, we describe how both the escape rate and low 
frequency current noise cause decoherence and report re- 
sults on measurements of these effects in Al/AlOx/Al 
Josephson junctions. 

Consider a Josephson junction shunted by ca- 
pacitance C, having a critical current To, and a parallel 
shunting impedance R(u>) due to the external wiring (see 
Fig. The supercurrent I through the junction is given 
by the Josephson relation I = Iq sin(7), and the voltage 
by V = ($o/27r)d7/dt, where 7 is the gauge-invariant 
phase difference between the superconducting wavefunc- 
tions on each side of the junction. For I < Iq, the phase 
7 may be trapped in a metastable well of the Josephson 
washboard potential, U = — ($o/27r)io COS7— (^o/2n)I^, 
or it may be in a running state with a non-zero average 
dc voltage 0. 

Quantizing the single junction system in the ab- 
sence of dissipation leads to metastable states that are 



localized in the wells (see Fig. EJ and adds the pos- 
sibility of escape to the continuum running states by 
quantum tunneling from the ith level with a rate T^oo. 
The energy barrier AU = (Iq^o/tt) ■ - [I/h) 2 ~ 
(7/7o)acos(///o)) to the continuum decreases as the bias 
current is increased, leading to a rapid increase in the 
tunneling rate with bias current |9j: 
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cillation frequency and iV s = AU/htOp is approximately 
the number of levels in the well. As the energy barrier is 
lowered, the energy of the states in the well move closer 
together and the well becomes more anharmonic until, at 
I = Io, the energy barrier disappears. 

The observed escape rate of the system from the 
zero-voltage state to the finite voltage state at a given 
bias point is T = UJLqIV-hx, Pi, where Pi is the probabil- 
ity of the junction being in the ith state. An ac current, 
I ac (either external or thermally generated) can induce 
transitions between levels i and j in the well with a rate 




FIG. 1: Circuit schematic of current- biased Josephson junc- 
tion. All the elements in the dashed box are represented by 
an equivalent resistance, R(lu). 
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K 1 2F^c(«'|7b')| ■ Since Ti^oo ~ 500^^00 for 
typical junction parameters, one expects to see a large 
enhancement in the escape rate if a microwave source 
is used to resonantly excite the system from the ground 
state |0) to the first excited state |1) (see Fig. EJ). 2] 

Each microwave resonance in this system will be 
broadened due to the interaction of the junction with 
noise transmitted via the wiring and described by the 
interaction Hamiltonian Hi nt — — (&o/2Tr)I no i se "f. Ther- 
mal noise and dissipation at the transition frequencies 
will cause changes in the populations of the states. At 
low frequencies, the resonant terms are insignificant and 
the noise only causes dephasing. 

Considering just the ground state |0) and the 
first excited state |1), transitions arise from thermal ex- 
citation from |0) to |1), a 1/RC decay rate from |1) to |0), 
and tunneling to the continuum, T^oo for i = and 1. 
At temperature T, the combined thermal and dissipative 
transition rates are:[10j 



RC(cxp(AE/kT) - 1) 



A_KJ RC(1 - exp(-AE/fcT)) v ' 

where AE — E\ — Eq is the difference in energy between 
the two levels. For kT <C AE, the upward thermal tran- 
sition rate is much smaller than the downward rate. From 
Eqn. n t ne tunneling to the continuum is much smaller 
for the ground state than the excited state in the an- 
harmonic region of interest where AU/huj ~ 3. 9] Thus 
we expect that the spectroscopic width of the |0) — > |1) 
transition is 

Ad! = Ti^oo + To^oo + To^i + Ti^o — Ti^o + Ti^oo 

~ 1/RC + ri^ 00 (4) 

Equations n and 01 imply that the level broadening, Auj, 
depends on bias through the ri_,oo term and should ex- 
ceed 1/RC as the bias current approaches Iq. 

To understand results on a real junction we must 
also take into account the dephasing effects of any cur- 
rent noise in the system. For sufficiently low frequencies, 




FIG. 2: Josephson junction potential energy U as a function 
of the phase difference 7. 



we can model this non-resonant decoherence mechanism 
as a simple smearing of the response with bias. This 
should result in a broadening of the spectroscopic width 
that depends on how sensitive the resonant frequency, to, 
is to changes in current, dui / dl . An rms current noise 07 
should produce an additional contribution to the spectro- 
scopic width of 2aidu}/dl . Including this current noise 
contribution in the previous form for the spectroscopic 
width gives: 

Au~l/RC + T 1 -> 00 +2(T I du/dI (5) 

Both the second and third terms in Eq. JSJ depend on 
bias, so that care must be taken in disentangling the two 
effects. 

Using double angle evaporation, we fabricated 
20 x 5 (urn) 2 Al/AlOx/Al Josephson junctions with 
J c ~ 14 A/cm 2 . Direct measurements of the junc- 
tion cur rent- volt age characteristics showed a subgap re- 
sistance of more than 10 4 VL at 20 mK. Escape rate mea- 
surements were made in an Oxford Instruments Model 
200 dilution refrigerator with a 20 mK base tempera- 
ture. We were able to tune the critical current of the 
junction by means of a superconducting magnet. The 
junctions were partially isolated from the bias leads by a 
10 nH surface mount series inductor and a 10 pF capac- 
itive shunt across the dissipative 50 f2 transmission line 
leads (see Fig. This isolation scheme was designed 
so that at the plasma frequency, the effective shunt resis- 
tance due to the leads would be stepped up from 50 f2 to 
much more than 10 3 Q, increasing the intrinsic Q of the 
system. To perform escape rate measurements, we start 
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FIG. 3: Enhancement of escape rate under 5.7 GHz microwave 
drive. Left axis is the difference in escape rate with and 
without microwaves divided by the escape rate without mi- 
crowaves. The large error bars on the left and right of figure 
come from a lack of counts in the escape histogram. The 
right peak is |0) — > |1) quantum transition, while left peak is 
|1) — » 1 2). Solid line is a Lorentzian fit to two peaks. 
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a timer and then ramp the current slowly (5 mA/s) using 
an HP 33 120 A function generator through a 47 kil resis- 
tor and monitor the junction voltage with a 2SK170 FET 
followed by an SRS560 amplifier. This output voltage is 
used to trigger the stop of timing, which is handled by a 
20 MHz clock. Escape events were binned in time with 
width t w ~ 50ns to create a histogram H(ti). The escape 
rate is then T(tj) = A- In [T,f =J H{U) /^^HiU)}. We 
convert the time axis to current by calibrating the ramp 
current as a function of time. 

We determine the spacing of the energy levels by 
comparing escape rate curves with (T m ) and without (To) 
a small microwave drive current applied. Figure [3] shows 
Ar/r = (r m — r )/r for a 5.5 GHz microwave signal. 
We chose the power so that Ar/r % 10 on resonance 
to ensure the occupancy of |1) is small. Two Lorentzian 
peaks are apparent, corresponding to the |0) — > |1) and 
1) — > |2) transitions. By measuring Ar/Fo for different 
applied microwave frequency we can measure how the 
bias current changes the energy level spacing of the |0) — > 
1} transition (see Fig. 

The data in Fig. also allows us to compute 
du/dl and convert the full width at half-maximum A J 
measured at each frequency (see Fig. 0Jd) to a width 
in frequency, Auj, or the spectroscopic coherence time 
associated with the two levels, r = 1/Aw. 

FigureElshows the coherence time r as a function 
of the center current of each |0) — * |1) peak. We note that 
the coherence time decreases markedly as J approaches 
Io ~ 14.12/iA, consistent with escape rate limiting of 
the lifetime of the upper state and excess low frequency 
current noise as in Eq. 5. 

In principle, it is possible for the effective shunt- 
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FIG. 4: (a) Drive frequency vs. center of the |0) — ► |1) reso- 
nance peak for Iq = 14.12/1 A The solid line is a smooth fit 
used only to extract a local slope. The dashed line is a fit to 
theory, (b) Full widths of each resonance for Jo = 14.12/iA. 
(Data set #050902) 



ing impedance R(w) to vary with frequency in such a 
way as to generate the changes in t(uj) seen in Fig. |S] 
We can rule out this explanation for the overall behavior 
of t(u>) by changing the critical current of the junction 
and remeasuring at the same frequency. Such a process 
changes T^oo but not R(u>) in Eq. 5. Results for two 
different Jo's are plotted in Fig. and^. Comparison 
of Figs. and 03 reveals that the coherence time at 
fixed frequency is lower for larger Jo. Since this measure- 
ment is at fixed frequency, the effect can not be due to 
R varying with frequency. 

To distinguish the effects of current noise and 
escape-rate broadening in Eq. we need to obtain an 
independent measure of the junction parameters. For the 
low critical current data, we fit the escape rate curves 
without microwaves 01 an d nn d Iq = 10.65 ± 0.01 fiA, 
C = 3.7 ± 0.3 pF, and T = 60 ± 3 mK. The 60 mK 
temperature was 40 mK above the base temperature, 
probably due to self-heating. We also numerically solved 
Schrodinger's equation (with hard wall boundary condi- 
tions) and chose Jo and C to fit the data in Fig. 
(dashed line). This yielded J = 10.645 ± 0.01 [iA and 
C = 3.7 ± 0.1 pF. The same analysis for the high Jo case 
gives J = 14.12 ± 0.01 [iA and C = 3.7 ± 0.1 pF. The 
quantitative disagreement in Fig. ^may come from not 
including corrections to the center peak locations due to 
current noise |l2| , the energy level shifting due to damp- 
ing, or a frequency dependent impedance (such as is sug- 
gested at 5.2 GHz in Fig. Eb). 

We now fit the coherence time data in Fig. 
by varying Jo and C and comparing the results with 
the previously determined parameters. We find Ti^oo 
by solving Schrodinger's equation numerically. To esti- 
mate the rms current noise, 07, we note that the full 
current width at half maximum shown in Fig. 0Jd never 
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FIG. 5: Coherence time r vs. bias current, I. Note that the 
escape rate from the ground state at 13.93 //A is around 10 3 /s 
while at 14.01/xA, it is around 3 x 10 6 /s. 
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FIG. 6: Coherence time, r, vs. bias current I. Solid lines 
are the theoretical fit for each data set. Lower frequency 
corresponds to larger current. The parameters for the fit in 
(a) are I = 14.12 fiA and C = 3.7 pF. For the fit in (b) 
la = 10.645 fiA and C = 3.7 pF. The dashed lines represent 
the contribution from the escape rate, and the dotted lines 
the contribution from current noise. 



drops below 10 nA. We thus assign 07 ~ 5 nA. To get 
a unique fit, we also assume RC 3> 1/ (2aidoj/dI). The 
solid lines in Fig. HJlshow the results of this procedure. 
The dashed lines show the contribution to the broaden- 
ing due to the escape rate alone, while the dotted lines 
represent the current noise contribution. The parame- 
ters for the lifetime fits, Iq = 14.12 fiA, C — 3.7 pF and 
Io = 10.645 fiA, C — 3.7 pF, agree with the parameters 
obtained from Fig. 01 verifying the inclusion of current 
noise and escape-rate-limited coherence in the model of 
Eqn. [SJ We note that as the bias current approaches Iq 
(low frequency) , the escape rate term begins to dominate 
the lifetime, while for lower currents (high frequency), the 
noise broadening dominates. 

To conclude, we have measured the resonance 
width of the transition between the lowest two quantum 
states in a Josephson junction qubit as a function of bias 
current, and found that the lifetime of the excited state 
falls rapidly as the bias current I approaches Iq . A model 
including continuous dephasing from tunneling as well 
as from current noise explains quantitatively the reduced 
coherence time. This ability to predict and calculate such 
junction behavior is crucial to the use of junctions in 



quantum computers and one of the reasons junctions are 
a good candidate qubit. 

For designs where low-frequency current noise is 
not a significant issue Q, consideration of the above re- 
sults in conjunction with Eqn. ^ suggests the following 
qubit design criterion. To obtain at least N op gate op- 
erations before decoherence occurs from tunneling alone, 
with each gate operation taking approximately N g ■ 2ir/to 
time, requires at least N s > ^ ln(N op N p ) + ^ ln(4327V s ) 



levels in the well. For N, 



op 



36 "n-"op-"p; 
= 10 6 and N n 



10, we find 



N s > 4. In the opposite limit, where current noise domi- 
nates, the junction must be biased at low currents during 
gate operations. 
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